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I. INTRODUCTION 



Lattice gauge theory has shown itself to be an extremely useful tool for studying 
non-perturbative physics. Its success is founded on the ability to systematically remove 
errors introduced by the discretization of space-time. The magnitude of such errors are 
determined by the lattice spacing a. In the continuum limit a these error terms vanish, 
and so one could expect to asymptotically approach continuum physics by moving to finer 
lattice spacings. Unfortunately the computational cost increases significantly as the spacing 
decreases. 



It is advantageous to 'improve' operators by identifying and algebraically eliminating 
discretization errors from their lattice definitions [|T|. In particular, calculations of the 
Yang-Mills action based upon the strategic combination of several different Wilson loop 
terms and the reduction of non-classical errors by tadpole improvement have demon- 
strated a significant decrease in deviations from expected continuum physics in smoothing 
algorithms such as cooling and smearing (see for instance Ref. ||^ and references therein). 
In this paper we consider the construction of (9(a^)-improved gluon field-strength tensors, 
F^y, from which both the topological charge density and the action may be constructed. In 
the following we refer to such an action as the "reconstructed action" . 



The accuracy of our improved operators is investigated by calculating the action, the 
reconstructed action, and the topological charge on an ensemble of field configurations. 
We then cool these rough configurations to produce highly self-dual configurations with a 
range of topological charges. Thermalization and cooling are carried out by use of parallel 
algorithms with appropriate link partitioning [^. Comparisons of the topological charge Q, 
and action normalized to the single-instanton action, S/Sq, provide quantitative measures 
of the merit of each of the improvement schemes considered. 



This paper is set out as follows: In Section II we look at the lattice definitions of the 
action and topological charge operators and give a brief discussion of the construction of 
the 0(a^)-improved forms of these operators In Section |T| we proceed to describe the 
construction of an analogously improved lattice version of the field-strength tensor, F^^. This 
tensor is used to construct an improved topological charge operator and the reconstructed 
action. In Section we use the cooling of gauge-fields on the lattice to determine which form 
of improvement gives the most continuum- like results, and compare the standard improved 
action with the reconstructed action to verify consistency. Our results and conclusions are 
presented in Section [V^. 



II. LATTICE ACTION AND TOPOLOGICAL CHARGE 



A. Wilson Action 



The lattice version of the Yang-Mills action was first proposed by Wilson . The action 
is calculated from the plaquette, a closed product of four link operators incorporating the 
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link U^, 

X fJ.<U ^ ' 

for an SU{N) field theory where the plaquette operator Wj^l^^^ is 

Wil""^^ = U^{x)U,{x + afi)Ul{x + av)Ul{x) 
The link variables U^lx) are defined by the Schwinger line integral 

Uf,{x) = VexpUg J dzAf,{z) \ , 



(1) 



(2) 



(3) 



and are in general non-Abelian (hence the V on the right-hand side of Eq. to denote 
path-ordering). Let us define Wj^^^'^^ to be the closed loop product in the fj, — u plane 
with extent m lattice spacings in the /^-direction and n lattice spacings in the i/-direction. 
Therefore Wj^J^^^ is the usual plaquette, W^]/^'^^ is an (a x 2a) loop 



^au"^^ = U^{x)U^{x + ajl)Uy{x + afi + ah')Ul{x + 2az>)f/J(x + ah')Ul{x) , 



(4) 



and so on. The Wilson action based on the plaquette for QCD (i.e., = 3) expanded 
around xq, the centre of the plaquette, may be written as follows 



'S'wil 



/5EE 

xo fJ,<U 

/?EE 

Xo fJ,<U 



1 



N 



ReTiWli'^'^Xo) 



1 - 1 + ^TrF^2,(xo) + Oig'a^') + 0{g'a') 



^TrF^2,(xo) + 0{g'a') + 0(^?V) 



EE + 0{a') + 0{g 



(setting P 



Xo ^<u 



EE 

Xo M:*^ 



(5) 



which reproduces the continuum action to O{o?). We note that the path-ordering is cru- 
cial to obtaining the non-abelian F^y = d^Ay—di,A^+ig[A^, A^] with errors of order 0{g'^a'^). 



In this investigation we use cooling to eliminate the high-frequency components of the 
gauge fields on the lattice. This reveals topological structures which correspond with classical 
minima of the action. The cooling algorithm works by calculating the local action at each 
lattice link, and updating the link to minimize this action. However the discretization errors 
inherent in the Wilson action cause the cooling algorithm to consistently underestimate the 
action at each lattice link 0. As a result the link updates do not accurately reflect the 
true structure of the gauge fields, leading to the destabilization of non-trivial topological 
structures (instantons and anti-instantons). 
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B. Improving the Action 



The standard Wilson plaquette action, S'wii, contains deviations from the continuum 
Yang-Mills action of order O(a^) and 0{g'^a'^). This problem may be remedied by improving 
the action. Tree-level improvement is a simple and effective means of eliminating lowest- 
order discretization errors by calculating the action from a combination of, for instance, the 
plaquette, 2a x a, and a x 2a rectangles. Since the plaquette and rectangles have different 
(9(a^) errors they may be added in a linear combination in such a way that the leading 
error terms cancel and one is left with the term corresponding to the Yang-Mills action plus 
(9(a^) and 0{g^a'^) errors. 

The form for a tree-level improved action can be easily determined from the quantities 

+0{g'a') + 0{g'a') , (6) 



li?eTriy(2xi)(^^) ^ l-^TrF2,(xo)-^TrF,.(xo)(492 + 9,2)F,.(a:o) 

+0{g\') + 0{g'a') , (7) 



\ReTTW^J<'\xo) = l-^TrF2,(xo)-^TrF,.(xo)(9j + 49^)F,.(xo) 

+0{g'a') + 0{g'a') . (8) 

Henceforth we shall make use of the definition 

P!ir''\x) = lReTTW^,T-\x). (9) 

We can therefore readily see that the 0{a'^) terms may be eliminated from the action by 
making the following construction 



"-^ X fl<U 
X 11,1/ 



TTF^^ixo) + Oia^) + Oig'a 



2Q V ^ -V 20 

(setting /3 = 



which reproduces the continuum action to 0{a ) and 0{g a 

This process eliminates classical error terms of order 0{a^) arising from the finite lattice 
spacing. Non-classical 0{g'^a^) errors arising from self-couplings of the gluon fields may be 
dealt with to some extent by the process of tadpole improvement 0. Due to the different 
numbers of links in plaquettes, rectangles, and other possible choices of Wilson loop, one 
redefines the value of each link by scaling it by the mean link uq (to account for the signifi- 
cant tadpole-style self-interactions of the gluon fields introduced on the lattice via the link 
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operators). In essence uq is an estimate of the sum of these unwanted tadpole terms, after 
the leading 0{g) terms in the links are factored out 0. Scaling the links by the mean link 

U,{x)^^, (10) 

Mo 

improves the accuracy of the expansions above. We use the plaquette measure 

which is updated after every sweep of cooling. When we include tadpole improvement 
factors, to correct for the difference between the numbers of links in the plaquette and the 
rectangular Wilson loops, the improved action takes the long-established form |^ 



5'lmp - o/^XH] 



(12) 

There is no reason why the tree-level improvement scheme needs to stop at this point. 
In principle it can be extended to arbitrary orders. For example, de Forcrand et al. |^ have 
used tree-level improvement to construct a lattice action which eliminates 0{a^) errors, by 
using combinations of up to five Wilson loop operators. These five loops, which we denote 
as Li, L5, have dimensions (1x1), (1x2), (2x2), (1x3), and (3x3). Of course, in the 
case of rectangular loops m 7^ n we average the contribution of the loops in each direction. 
Hence we define 





- p(lxl) 

— ? 




L2 


— p(2x2) 




L3 


= l/p(2xl) 
— 2 I M*^ 


1 p(lx2) 




= 1 fp(3xl) 
— 2 I 


1 p(lx3) 




— p(3x3) 





The improved action of de Forcrand et al. can then be written in the form 

^ 1 

^ = E7-ir^c,5,, (13) 



i=l [^n- 



where Si is the action calculated by inserting the loop Li in place of Uf^u in the definition 
of the Wilson action, Eq. (|I|), where m and n are the dimensions of the loop Li, and where 
the Ci are improvement constants which control how much each loop contributes to the total 
improved action. In our work we explicitly include tadpole improvement factors, 

(14) 



(2m+2n-4) ' 
Mo 



whereas de Forcrand et al. left these factors as unity in Eq. (pIBf). Our experience is 



that tadpole improvement factors are beneficial in the early stages of cooling and remain 
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beneficial even as mq — 1. 



The values of the improvement constants which eliminate order 0{a^) tree-level error 
terms are 

Cl = (19-55c5)/9, 
C2 = (l-64c5)/9, 
C3 = (640c5 - 64)/45, 
C4 = 1/5 -2c5, 

where C5 is a free variable. By adjusting C5 we create different 0(a'^)-improved actions, 
which in general will have different 0{a^) error terms. This means that we can attempt to 
minimize 0{a^) errors by tuning C5. If we set C5 = 1/10, then C3 = C4 = creating a "3-loop" 
action, i.e., an action constructed from three of the five loops. Setting C5 = we create a 
"4-loop" improved action. To construct a specific 5-loop action, de Forcrand et al. chose 
to set C5 = 1/20, halfway between the 3-loop and 4-loop values to stabilize single instanton 
sizes under coohng. These so-called 3-loop, 4-loop, and 5-loop improved actions are all free 
of 0{a^) errors, and have been shown to preserve (anti-)instantons for many thousands of 
sweeps. 



C. Topological Charge 

On the lattice, the total topological charge is obtained by summing the charge density 

^^^^ ^ ^^e^^P<xTr{F^^(a;)Fp^(a;)} , (15) 

over each lattice site 

Q = E?(^) > (16) 

X 

where /i, u, p, a sum over the directions of the lattice axes. Further to their construction 
of the improved action, de Forcrand et al. have used an analogous procedure to define an 
improved lattice topological charge operator |^. We choose to construct the C(a'^)-improved 
topological charge in a different manner. In the next section we derive an improved version 
of the field-strength tensor, F^^. This improved operator forms the basis of our definition of 
the improved topological charge. The derivation of the improved field-strength tensor will 
serve to illustrate the general methods used to deduce the results already stated (without 
proof) above for the improved action. 



III. THE LATTICE FIELD-STRENGTH TENSOR 



Consider the problem of expanding a generic Wilson loop operator. Such an operator is 
defined as the path-ordered exponential of the closed path integral around the loop 

p|/{mxn) _ p ^ig§A.dx^ ^yj^^ 
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As noted earlier, the path ordering is crucial to obtaining the full non-abelian field-strength 
tensor in the action associated with W^^^") such that the errors are 0{g^a'^). Note that 
mean-field improvement reduces the magnitude of these O{g'^o?) errors present early in the 
cooling process. The 0{g'^a?) errors are also suppressed via cooling and we do not consider 
them further. 

We are considering tree-level improvement to remove 0{a^) and higher errors. It is 
therefore sufficient to work to leading non-vanishing order in g, which is equivalent to the 
abelian theory. We hence determine the coefficients in the expansion of any Wilson loop by 
considering the weak coupling limit g ^ 

Rewriting the abelian version of the line integral / A.dx in Eq. (0) as §{A^dx^ + Aydx^,) 
we see that we can apply the two-dimensional version of Stoke's theorem (i.e., Green's 
theorem) 

r r r ( dK dJ\ 

f Jdx + Kdy = / — — dxdy , (18) 

JdR J J R \ OX oy I 

where OR is the closed contour containing the two-dimensional surface R. Therefore we see 
that we can readily rewrite our integral around the loop as 



//•ma/2 /•na/2 
A.dx = / dXf, dx^ {d^A„{xo + x) - d^Af,{xo + x)) . (19) 

J— ma/2 J— na/2 



where xo is the centre of the loop. The integrand is identified as the abelian field-strength 
tensor 

//•ma/2 /•na/2 
A.dx = / dx^ dx^F^,y{xo + x) . (20) 

J— ma/2 J— na/2 

Hence the expansion of any abelian loop will be in powers of F^^,, reproducing the F^j^F^^'^ 
term of the action at order Olg"^), and containing error terms of relative order 0{g'^a'^) and 
O(a^). 

Now we Taylor expand the F^^^s around the point Xq 

//•ma/2 /•na/2 I 

A.dx = / dx^ dx,y iFf,^{xo) + ^XadaFf,^{xo) 

J— ma/2 J— na/2 \ q, 

+ ^Z!^"^/3^«^/3^m;^(^o) + ■ ■■ I , (21) 

Since our plaquette is a two-dimensional object we sum over the /x and z/ directions, {a = 
or z/, and similarly for j3). 



//•ma/Z rna/Z / 

A.dx = / dXf, dx,y i F^,y{xo) + {Xf,df, + x,ydu} Ff,^{xo) 

J— ma/2 J— na/2 \ 

+ {x^x^df^d^} F^^{xo) + - {xldl + xldl + d^dl^F^^,{xQ) 
+ \ + + ^l^-^l^- + F^.(xo) 



+ \ {^Ixldldl + xlxldldl] F,^{xo) +0{x')) . (22) 



The terms with odd powers of a;^ or x^, vanish upon integration, because of the symmetric 
integration hmits, about the center of the loop. Let us perform the expansion of the standard 
(1 X 1) plaquette. 



^ A.dx = j j dx^dx^ (^F^u{xo) + ^ + F^^{xo) 



= a^F^^ixo) + ^{dl + dl) F^.(xo) 

+ {dl + dt) F,^{xo) + ^ {dpi) F,Axo) + 0{a') . (23) 

This expansion can be determined to arbitrary order by taking the Taylor expansion in 
Eq. ( piD to the desired order. 

So consider again the form taken by the equation for the expansion of the (1x1) loop 
operator in the non-abelian case, which derives from the Schwinger line integral around a 
closed path 



W;(ixi) = 'Pe'^^A-'i- = p |i + A.dx - ^^{j>A.dxf + ■ ■ - j . 



(24) 



The standard approach to constructing the action is to take the real part of the trace of this 
expansion, thereby extracting the leading term which is just one, and the third term which 
is proportional to F^^^. However we are interested here in extracting the term proportional 
to Ffj,y. We can do this by making the following construction, which takes advantage of the 
hermitian properties of the Gell-Mann matrices 



1x1) = r^l + igjAdx-^^{jAdxf + 0{g')^ , 
iy(ixi)t = V^l-igjAdx-^^{jAdxf + 0{g')^ , 



and hence 



—I 

T 



■^(1x1) _ ^(ixi)t _ lTr(H^(ixi) - W^(ixi)t)l ^ ^Adx + 0{g 



= ga^ [f^.(xo) + 0{a^) + 0{g^a^) 



(25) 



We have subtracted one-third of the trace to enforce the traceless property of the Gell-Mann 
matrices. 
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Eq. ( p5D demonstrates how we construct the field-strength tensor F^i, from the (1x1) pla- 
quette. Since the expansion in Eq. (^) corresponds with a given Wilson loop only through 
the choice of integration limits, this suggests that it is also possible to construct an improved 
field-strength tensor from a suitably chosen linear combination of Wilson loops. This im- 
proved field-strength tensor can be fed directly into Eq. (|T3|) resulting in a topological charge 
which is automatically free of discretization errors to the same order in a^. Furthermore, 
since the Yang-Mills action is based upon the field-strength tensor (squared), it is possible 
to create a reconstructed action, Sr, by inserting the improved field-strength tensor directly 
into the equation 

^i^ = /?EE^TrF,.(x)^ (26) 

X fJ,,U 

Since the reconstructed action is improved differently to the standard improved action 
operator, Eq. ([I3|), and constructed from clover terms, the value of the action calculated 
with the reconstructed action operator can be compared with the value calculated with the 
standard improved action operator at each cooling sweep as a mechanism to explore the 
size of the remaining discretization errors. 



IV. IMPROVING THE FIELD-STRENGTH 



We now wish to construct sums of clover terms designed to remove O(a^) and 0{a^) 

jjLU • 



errors relative to the leading term, a?Fnu- Let us begin by defining 



A = ga^F. 



l^U 1 

B = ga' {dl + dl) F,, , 

flU 1 



C = ga' + dt) F, 

= 9a' {dldl) F,^ . (27) 

Let us now denote as the combination of Wilson loop terms extracted from the 

left-hand-side of Eq. ( pSj) corresponding with the loops used to construct a clover term as 
depicted in Fig. |I]. Constructing the four clover terms and symmetrizing in n ^ m we have 

If rma rna rO rna 

C(™'") = ^ \ dxj dx,F^, + / dxj dx,F^, 

10 Jo J-ma Jo 



+ 



/O rO pma rO 

dx ^ I dXi/F^jy -\- I dx ^ I dXyF^if 
-ma J —na Jo J —na 

rna rma rO rma 

~\~ I dXn j dXyFuy "j" / dX n I dXyFui/ 

Jo Jo J-na Jo 



/O rO rna rO 

dx^ / dXyF^^ + dXf, dx^F^ 
-na J —ma Jo J —ma 

r rma rna rna rma '\ 

— S / dX^ / dXjyF^i/ -|- / dX^ / dXyF^y (28) 

O vJ —ma J —na J —na J —ma ) 



Expanding the loop operators as in Eq. (|2^ ) we find that 

^ ^+6^120^36 
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/ 



8 



->- 



/\\/ 



->- 



->- 



->- 



FIG. 1: The four m x n loops and four n x m loops used to construct the clover term from which 
we construct the lattice topological charge operator. 



(7(2,2) 






16^ 


C(l,2) 


D 






(7(1,3) 


= 3A+^B + 






(7(3,3) 


27 

= 9A+—B^ 


40 





Notice that the C(a^) terms in each F^i, will be "promoted" to C(a^) (and higher) by 
taking the product F^i,Fp„, while the 0{a^) terms will be promoted to 0{a}^) (and higher). 
Consequently, since we wish to eliminate 0{a^) errors from the product, we need only 
expand out to (and eliminate) order O(a^) corrections in F^^. We therefore have five 
equations and four unknowns. 



The improved field-strength may be written 

i^ir - hC^^^ + hC^f^ + fcaCiL'^) + hCi^'^ + hC^f , (30) 

where the ki, the "weightings" of each loop, are defined by 

hA^^'^^ + A;2>1(2'2) + hA^^^''^ + k^A^^^^^ + A;5^(=^'3) 
A;iB(i'i) + A;2H(2.2) + k^^B^^,-^) + k^B^^^"") + hB^^^^^ 

k^C^^'^) + A;2C(2.2) + A;3C(i'2) + k^C^^^^) + fcsC^^-s) 

such that the weighted coefficient of ga^F^i, is 1 and the weighted coefficients of the 0{a?) 
and 0{a^) terms vanish in the sum of loops. Here ^(™'") represents the coefficient of the 
A term from the expansion of C^^'"^^ and so forth, with the values of the coefficients taken 



1, 
0, 
0, 
0, 



(31) 
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from Eqs. p9|). To find the values of the improvement constants we use these coefficients to 
construct an equivalent matrix equation. Eq. (pl|) takes the form 



1 


4 


2 


3 


9 


1 


8 


5 


5 


27 


6 


3 


6 


2 


2 


1 


8 


17 


41 


243 


120 


15 


120 


40 


40 


1 


16 


2 


3 


81 


36 


9 


9 


4 


4 



ki 



Using Gauss- Jordan elimination, the improvement coefficients are 



ki 


= 19/9- 


55 /cs , 


k2 


= 1/36- 


I6/C5 , 


k3 


= 64A;5- 


32/45 , 


^4 


= 1/15- 


6/C5, 



(32) 



where the coefficient of the 3x3 loop (in this case k^) is a tunable free parameter. 



We can see that if we set k^ = 1/90, then k^ = k^ = 0, eliminating the contribution from 
the C^^'^-* and C^^'^-* loops and creating a 3-loop (9(a^)-improved field-strength tensor. We 
may create a 4-loop improved field-strength tensor in three different ways, by setting fcs = 0, 
/cs = 19/495, or k^ = 1/576, where the 3x3, 1x1 or 2x2 loops are eliminated respectively. 
We focus on a 4-loop improved tensor with ^5 = throughout this investigation, analogous 
to the 4-loop action of de Forcrand et al. 



It must, of course, be noted that the 3-loop and 4-loop action and field-strength operators 
are simply special cases of the 5-loop operators, corresponding with particular choices of the 
parameters C5 and k^ respectively for the action and field-strength tensor. As noted above, 
we have already chosen to investigate a 5-loop action by using de Forcrand et al.'s choice of 
C5 = 1/20, halfway between the 3-loop and 4-loop values of C5. In an analogous manner we 
select the value k^ = 1/180 (midway between the 3-loop and 4-loop values) for the 5-loop 
improved field-strength tensor. 



V. COMPARISON OF IMPROVEMENT SCHEMES 



The configurations used in this work are constructed using the Cabbibo-Marinari p| 
pseudo-heatbath algorithm with three diagonal SU{2) subgroups looped over twice. We 
thermalize for 5000 sweeps with an (9(a^)-improved action with tadpole improvement from a 
cold start (all links set to the identity) and select configurations every 500 sweeps thereafter. 
Cooling proceeds via updates of the three diagonal SU{2) subgroups, looped over twice 0. 
Configurations are numbered consecutively in the order that they were saved during the 
process of thermalization. Hence configuration 1 was saved after 5000 thermalization sweeps 
from a cold start, configuration 2 was saved 500 sweeps after configuration 1, configuration 
3 was saved 500 sweeps after configuration 2, and so on. Our results are generated on a 
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12^ X 24 (untwisted) periodic lattice at (3 = 4.60, with a lattice spacing of a = 0.125 fm. 



As a thermalized configuration is cooled over the course of many hundreds of sweeps, 
the action of the configuration monotonically decreases. This occurs because the cooling 
algorithm smooths out short-range fluctuations in the field. As these ultraviolet components 
of the field are suppressed, the underlying medium and long-distance structure of the field 
is revealed. When there are no regions of positive and negative topological charge density 
in the process of annihilating, the configuration can be thought of as becoming a dilute 
instanton gas. If cooling proceeds for long enough, the only non-trivial parts of the field 
which will remain are instantons or anti-instantons. 

The total action and topological charge of a configuration satisfy the condition ||^ 

S>\Q\So = \Q\^. (33) 
9 

Since all continuum single (anti-)instantons are non-trivial minima of the local action with 
\Q\ = 1, they saturate the equality in Eq. ( p3D and therefore have an associated action of 
5*0 = Svr^/^f^. For well-cooled configurations it follows that we may estimate the number of 
instantons on the lattice by dividing the total action by the quantity 5*0 = Sn'^/g'^. We shall 
see that one of the principal criteria by which we may judge the value of an improvement 
scheme is how closely the values of 5*/ 5*0 and the topological charge approach integer values 
as we continue cooling and approach self-duality. 



In a sufficiently dilute instanton gas the only contributions to the action and topological 
charge of the field will be that due to the (anti-) instantons. Hence if there are nj instantons 
and ua anti-instantons well-separated on a large volume, the total action and topological 
charge should be integers, and satisfy the conditions S/Sq = rij + ua and Q = rii — ua- As 
cooling proceeds, instanton- anti- instanton pairs will annihilate until S/Sq = \Q\ at which 
point the configuration is self-dual. What we expect to see is that as we begin cooling, Q 
should relatively quickly become integer and remain stable thereafter. As we continue to 
cool towards self-duality we should then see S/Sq monotonically decreasing to the point 
where S/Sq — > \Q\ from above. 



Our lattice volume is not large enough that we observe configurations, stable under 
cooling, with both instantons and anti-instantons simultaneously, i.e., on our lattice they 
annihilate as cooling proceeds. In Fig. ^ we present the action for two configurations 
each having a topological charge (stabilized after preliminary cooling) of |Q| = 1. As 
cooling proceeds one observes that S/Sq repeatedly decreases by an increment of two. 
During this period the topological charge remains stable. This indicates that the action 
is dropping due to instanton-anti-instanton annihilation events. Notice that the plateaus 
persist for longer as cooling continues. It is known that there can be no fully self-dual 



\Q\ = 1 configurations on the torus JTO], [Tl|]- The fate of = 1 configurations under 



continued cooling will be discussed elsewhere [0. We have observed that configurations 
which stabilize at larger topological charges do not exhibit such pronounced plateauing. 
It is reasonable to deduce that as the lattice becomes more sparsely populated, the mean 
period between instanton-anti-instanton annihilation events becomes longer, as it becomes 
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500 1000 1500 2000 2500 

Sweeps 



FIG. 2: Action values of configuration 1 (thin line) and configuration 11 (bold line), cooled with 
the 3- loop improved action. Both configurations have a total topological charge of 1. It can be 
clearly seen that the action of each configuration repeatedly plateaus briefly near integer values 
(dashed lines) before dropping by two units. This is characteristic of instanton-anti-instanton 
annihilation events. In the case of configuration 11, plateauing begins around sweep 500 {S/Sq ~ 9). 

possible to have well-separated, weakly-interacting instantons and anti-instantons on the 
lattice. On a sufficiently large lattice, for sufficiently small total \Q\ it should then be 
possible to arbitrarily closely approach a true dilute instanton gas with both instantons and 
anti-instantons, where the configuration is locally (but not globally) self-dual. 

On the lattice, discretization errors are expected to lead to non-integer values of S/Sq 
and topological charge for a highly-cooled configuration. The better the improvement, the 
closer these values will come to integers. Fig. ^ shows how S/ Sq and the topological charge 
approach the same value, indicating the approach to a self-dual configuration. They remain 
stable for several hundred sweeps of improved cooling. We now wish to determine which 
cooling scheme (3-loop, 4-loop, or 5-loop) and which field-strength improvement scheme are 
best. 

Throughout the discussion of our results we shall use the following notation: the cooling 
action will be denoted by S, the reconstructed action by Sr and our highly-improved topo- 
logical charge by Q. In some cases the type of improvement scheme used will be denoted 
by a number in parentheses. Hence the 3-loop improved cooling action is written as 5'(3), 
our topological charge calculated from a 4-loop improved field-strength tensor is written as 
Q(4), our 5-loop reconstructed action is written as Sr{5), and so on. 
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FIG. 3: An example plot of how improved cooling stabilizes the action (circles) and topological 
charge (triangles) at consistent values. The action is rescaled by dividing by Sq, the action of a 
single instanton. Cooling is performed with an S'(3) action, and the topological charge is (5(3). 
[The notation used is that S{n), SR{n), and Q{n) are the n-loop forms of the improved action, 
our reconstructed action, and our improved topological charge operator respectively.] 

A. Improved Cooling 

We shall commence the assessment of which improved action has the smallest discretiza- 
tion errors by demonstrating explicitly the discretization errors in the standard Wilson 
(plaquette) action. Fig. ^ shows the action against sweep number for a configuration 
cooled with the Wilson action (triangles), denoted as 5(1), and the S{2) action (circles) 
from Eq. (p!2D. It can be clearly seen that the Wilson action drops to a temporary 
plateau, but eventually destabilizes and drops by approximately one unit around sweep 
number 250, and again around sweep number 750. It should also be noted that when the 
numerical data are examined closely these plateaus are found to occur somewhat below 
integer values. However the S{2) cooling scheme plateaus at a value much closer to inte- 
ger (in this case 6) and remains at this value without destabilizing after 1000 cooling sweeps. 

This result (and other similar results on other configurations) are an indication that 
the ordinary Wilson plaquette action significantly underestimates the true action of the 
configuration. Because the large discretization errors in the action prevent the cooling 
algorithm from updating the links in a manner which accurately reflects the classical 
structure of the fields, plaquette cooling destroys important topological structures. It 
is true that the action may be expected to drop (while the topological charge remains 
constant) due to instanton-anti-instanton annihilation during the cooling process. However, 
since the action drops by an increment of one unit (not two) it is clear that we are not 
observing a decrease in the action as a result of instanton-anti-instanton annihilation, but 
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FIG. 4: Action of configuration 89 cooled with S{1) (triangles) and S{2) (circles) cooling schemes 
after 1000 sweeps. [See Fig. (^) caption for notation.] 



rather the destabihzation and destruction of a single instanton or anti-instanton. 

While the S'(2) results are promising, it has already been shown analytically that 5'(2) 
does not stabilize instantons on the lattice 0. Removal of O(a^) errors and long-term 
stabilization of instantons requires the consideration of additional loops. 

We now shift our attention to comparisons between the various improved actions. Figs. |^ 
and P show the value of the action attained by two different initial hot configurations as 
they are cooled with various improved actions. We can see that in both cases the 5'(2) 
action plateaus well below the relevant integer value, due to large 0{a'^) discretization 
errors with negative sign. We furthermore notice that in the case of Fig. ^, the >S'(4) action 
drops below the integer plateau at S/Sq = 6. In fact, while not shown here it continues 
falling to a plateau at S/Sq = 5, indicating that it cools away an instanton-like structure 
which the other actions do not. When we examine Fig. ^ we see that the discretization 
errors in the 5(4) action are relatively large and have negative sign. This suggests that 
the 5(4) action consistently underestimates the action of the configuration to which it is 
applied. In the case of Fig. ^ this underestimation was severe enough to completely destroy 
an instanton during the cooling process. This in itself is quite a remarkable result, since it 
demonstrates that 0{a^) discretization errors can be large enough to destabilize instantons 
in some cases. We shall henceforth discard 5(4) cooling on the basis that it has large 0{a^) 
discretization errors with negative sign. 
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B. Improved Topological Charge and Reconstructed Action 

In order to determine which improvement scheme for F^^, produces the most continuum- 
like results we will assess which topological charge operator produces results closest to 
integer values. The top diagram of Fig. |^ shows the development of the topological charge 
of a hot configuration over the course of 2400 cooling sweeps with the 5'(3) improved 
cooling action. The Q(l) topological charge is too far from an integer value to be seen on 
the scale we have chosen for this diagram, but we can clearly see that the (5(3) and Q(4) 
operators produce significantly better results than the Q{2) operator. At the bottom we 
see equivalent results for configuration 89 (used in Fig. over 1200 sweeps. This time we 
see that the (5(3) operator gives marginally more integer-like results than those observed 
with the (5(4) operator. But again the QiX) and Q{2) operators are substantially more 
inaccurate, as we would expect from the order of improvement used to construct (5(3) and 
(5(4). Our (5(5) operator, constructed with a value of midway between the 3-loop and 
4-loop values, produces near-perfect results for both parts of Fig. 0, sitting between the 
(5(3) and (5(4) values. Further fine tuning of appears unnecessary as (5(5) lies above and 
below the integer values in the top and bottom plots of Fig. ^ respectively. 

Since the improved field-strength tensor plays no role in affecting the structure of the 
field configurations as they cool, we consider an ensemble of configurations that have 
already been cooled to self-duality (for several hundred sweeps) and assess their topological 
charge (and reconstructed action) with each different improved F^^. Table | shows these 
results for a number of configurations, all >S'(3) cooled. We can see that in three of the 
four cases the 5-loop improved topological charge, (5(5), and reconstructed action, 5^(5), 
produce results that are both closest to integer values, and in closest agreement to the value 
of the cooling action. In the case of configuration 56, Sr{?>) and (5(3) produce the results 
closest to integer values. 

A further comparison, between results for 5'(3) and 5(5) cooling, is presented in 
Table |T|. In this case we wish to determine whether the choice of cooling action affects the 
dependability of the reconstructed action. We have not chosen to consider the 5(4) cooling 
as it has already been deemed unreliable. We see that in the later stages of cooling, S^ib) 
gives values closer to the cooling action used in each case (whether the cooling action is 
3-loop or 5-loop), which suggests that it may be a more accurate probe of the structure of 
the fields produced by the cooling algorithm than either the 3-loop or 4-loop reconstructed 
action operators. For these reasons it appears that the 5-loop improved field-strength 
tensor with k^ = 1/180 is the most dependable of the choices we have examined. 



VI. DISCUSSION 

We have now considered the relative merits of the various improvement schemes for the 
cooling action. We have also considered the relative merits of the various improvement 
schemes for the reconstructed action and the topological charge operators (which of course 
are indicative of the accuracy of the improved field-strength tensor). It appears that 
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FIG. 5: Action values over the first 1000 cooling sweeps of configuration 89 cooled with the 5(2) 
(crosses), 5(3) (triangles), 5(4) (vertical line) and 5(5) (circle) actions. Note that 5(4) drops to a 
value near 5.00, while the other actions plateau near 6.00. [See Fig. (H) caption for notation.] 
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FIG. 6: Action values over the first 2600 cooling sweeps of configuration 32 cooled with (from top 
to bottom on the right-hand-side of the figure) 5(3) (triangles), 5(5) (circles), 5(4) (squares), and 
5(2) (crosses). [See Fig. (^) caption for notation.] 
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FIG. 7: Topological Charge of configuration 27 (top) and configuration 89 (bottom) cooled exclu- 
sively with the S{3) cooling scheme, for 2400 sweeps and 1200 sweeps respectively. In descending 
order the curves are Q(3), Q(5), Q(4), and Q{2). The Q(5) (circles) operator falls directly between 
the Q{3) (triangles) and (5(4) (squares) topological charge operators, at almost perfect integer 
values. [See Fig. (0) caption for notation.] 
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TABLE I: Reconstructed action (rescaled by ^o) and Q values from various configurations with >S'(3) 
cooling. The columns are the configuration number, the sweep number at which the calculation 
was performed, the value of the cooling action, the values of 5/j(2) to Sr{5), and the values of 
(5(2) to (5(5). [See Fig. (^) caption for notation.] 



Config. 


Sweep 


S{3) 


Sr{2) 


Sr{3) 


Sb.{4) 


SRi5) 


89 


1071 


5.99992 


5.98908 


6.00110 


5.99777 


5.99943 


27 


1400 


2.00013 


1.99687 


2.00055 


1.99965 


2.00010 


90 


1046 


1.99999 


1.99849 


2.00012 


1.99978 


1.99995 


56 


2446 


1.99996 


1.99878 


2.00002 


1.99979 


1.99990 



Config. 


Sweep 


5(3) 


Q(2) 


Q(3) 


Q(4) 


Q(5) 


89 


1071 


5.99992 


5.98899 


6.00109 


5.99777 


5.99943 


27 


1400 


2.00013 


1.99680 


2.00050 


1.99961 


2.00006 


90 


1046 


1.99999 


1.99848 


2.00011 


1.99977 


1.99994 


56 


2446 


1.99996 


1.99878 


2.00002 


1.99979 


1.99990 



the 5-loop improved field-strength is most accurate, while the 3-loop field-strength is 
substantially cheaper (in a computational sense) and still produces excellent results. 



Assessment of improvement in the cooling action is more difficult. It has been stated 
by de Forcrand et al. that S{^) cooling is most dependable and produces the best 
stabilization of instantons. We have cooled configurations with a topological charge of 
|(5| = 2 and have seen that the instantons rapidly settle on mutually consistent sizes, and 
furthermore each instanton changes size by a factor of no more than 1.1 over the course 



of several thousand sweeps. This work will be reported in detail elsewhere |T^. However 
de Forcrand et a/.'s work establishes a precedent which must be recognized, so let us 
mention some differences between their work and ours. Firstly, our work is performed on 
configurations which contain several instantons or anti- instantons while de Forcrand et a/.'s 
work is performed on configurations which consist of a single instanton, obtained via 5'(1) 
cooling. Furthermore our work is performed on a 4-toroidal mesh of lattice points, with 
untwisted boundary conditions while de Forcrand et al.^s work is performed with twisted 
boundary conditions. The twist in the boundary conditions serves to stabilize the single 
instanton present, since self-dual |(5| = 1 configurations are not permitted on an untwisted 



torus TO, 11, 13 



As demonstrated in Table ||, for our (periodic) lattice the S{3) cooling action produces 
results as good as S{5). S{3) cooling is significantly faster than S{5), since it requires fewer 
link multiplications, and gives no noticeable decrease in long-term instanton stability. In 
fact, for large numbers of cooling sweeps S{5) cooling appears to consistently produce lower 
values of S/ Sq than S{3) cooling, and lower than the required integer value, suggesting the 
presence of 0{a^) discretization errors which contribute negatively. Over the long term 
these errors may destabilize instantons (as was observed with S{4) cooling), suggesting that 
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TABLE II: Reconstructed action and cooling action values (all rescaled by So) from various con- 
figurations cooled with S{3) (columns 3 to 5) and S{5) (columns 6 to 8). [See Fig. (^) caption for 
notation.] After no more than 170 cooling sweeps, the mean-field improvement factor uq for con- 
figurations 89 and 32 have settled on a value of 0.99992, and remain completely stable thereafter, 
while configuration 39 settles on a value of mq = 0.99996. The resultant weighting of the 3x3 loop 
to the 1 X 1 loop is 0.9995. 



Config. 


Sweep 


5(3) 


Sr{3) 


Sr{5) 


5(5) 


Sr{3) 


Sr{5) 


89 


200 


6.03297 


6.03417 


6.03243 


6.05030 


6.05190 


6.05012 




600 


6.00010 


6.00130 


5.99961 


5.99957 


6.00116 


5.99945 




1000 


5.99992 


6.00110 


5.99944 


5.99956 


6.00115 


5.99945 


39 


220 


3.12225 


3.12234 


3.12203 


3.16616 


3.16629 


3.16596 




520 


3.00024 


3.00041 


3.00013 


3.00051 


3.00074 


3.00045 




920 


2.99992 


3.00009 


2.99980 


2.99986 


3.00009 


2.99980 


32 


500 


6.00061 


6.00300 


5.99995 


6.00000 


6.00319 


6.00006 




1000 


6.00051 


6.00284 


5.99988 


5.99984 


6.00298 


5.99990 




1500 


6.00047 


6.00274 


5.99986 


5.99983 


6.00294 


5.99989 



S{3) is the safer choice for investigations of long-term instanton stability. 



VII. CONCLUSIONS 



In summary, to assess the relative merits of the improvement schemes discussed in this 
paper, we compare the various cooling actions (rescaled by the single instanton action 
5*0), with reconstructed actions (rescaled by So) and topological charges based upon 
various forms of the improved field-strength tensor, F^^,. For studies of instantons on 
periodic (untwisted) lattices the S{3) action is optimal for cooling, due to the superior 
accuracy with which it approaches integer values, its promise of long-term stability, and 
cost-effective execution. Five-loop improvement is optimal for the field-strength tensor, 
as most integer-like results are obtained from Q{5) for the improved topological charge 
operator and Sr{5) for the reconstructed action. However the 3-loop operator also produces 
excellent results, and is much less computationally demanding than the 5-loop operator. 
This recommends the 3-loop field-strength as a reasonable alternative choice to the 5-loop 
improved field-strength when computational resources are limited and improved accuracy 
at the order of a few parts in 10^ is not required. 



The improved field-strength tensor has already been used in a number of investigations. 
For example, the 3-loop topological charge described in this paper has been used to assess 
the Atiyah-Singer index theorem \H\ on uncooled and cooled field configurations [^. The 
highly-improved field-strength tensor has also been used in a "fat-link irrelevant clover" 
(FLIC) action |]l6l, a modification of the Sheikholeslami-Wohlert improved quark action 
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T7| , where the hnks of irrelevant operators are fattened via APE smearing. Detailed 
investigations of the long-term behaviour of instantons under (9(a^)-improved cooling have 
been performed [|12], including investigations of the consequences of the Nahm transform 
||T0| , |T3|, |18| for the stability of |(5| = 1 configurations on the 4-torus. 

Areas of future research using the improved field-strength tensor include adaptation to 
the study of glueball and hybrid meson spectra where colour-magnetic and electric fields 



explicitly provide gluonic excitations in the hadron interpolating fields |]19|, Future 
investigations should also investigate dislocation thresholds to learn how coarse the lattice 
can be made before a substantial reduction in the quality of data occurs. 
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